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Most of the baryons in galaxy clusters reside between the galaxies in a hot, 
tenuous gas 1 . The densest gas in their centers should cool and accrete onto 
giant central galaxies at rates of 10 — 1000 solar masses per year 1 . However, no 
viable repository for this gas has been found 1 . This paradigm changed abruptly 
when new X-ray observations showed far less cooling below X-ray temperatures 
than expected 2 . Consequently, most of the gas must be heated and maintained 
above ~ 2 keV 3 . The most promising heating mechanism concerns powerful 
radio jets emanating from supermassive black holes in central cluster galaxies 4 . 
Here we report the discovery of giant cavities and shock fronts in a distant 
(z = 0.22) cluster caused by an interaction between a radio source and the hot 
gas surrounding it. The energy involved is ~ 6 x 10 61 erg, the most powerful radio 
outburst known. This is enough energy to quench a cooling flow for several Gyr, 
and to provide ~ 1/3 keV per particle of heat to the surrounding cluster. 

Cavities with diameters ranging from a few to a few tens of kpc have been found in the 
hot gas surrounding nearly two dozen galaxies, groups, and clusters 5 . Their enthalpy (free 
energy), which ranges between pV = 10 55 — 10 60 erg, scales in proportion to the cooling X-ray 
luminosity and the radio power of the host system 5 . The cavities and weak shocks in half 
of these systems are currently injecting enough energy into the hot gas to balance radiation 



losses (cooling) 6 ' 7 , but it is unclear whether they can quench cooling over longer timescales. 
Systems without cavities today could have been heated in the past by powerful but relatively 
rare outbursts 8 . However, prior to the discoveries of large scale cavities and shocks discussed 
here and in a companion paper 9 , little evidence existed to support this conjecture. 

Giant cavities, each roughly 200 kpc in diameter, were found in a Chandra X-ray image 
the optically poor cluster MS0735. 6+7421 (Fig. 1). The center of the cluster harbors a cD 
galaxy that hosts a radio source roughly 550 kpc in size. The radio lobes fill the cavities 
suggesting, as in other clusters, that the gas was displaced and compressed by the advancing 
radio source. Both the cavities and the radio source dwarf the cD galaxy, which is itself a 
member of the class of the largest galaxies in the Universe. 

The average pressure surrounding the cavities (Fig. 2) is P ~ 6 x 10~ n erg cm~ 3 . The 
work required to inflate each cavity against this pressure is pV ~ 10 61 erg. Their enthalpy can 
approach ApV per cavity, depending on the equation of state of the gas filling them, giving 
a value of 8 x 10 61 erg. To place this figure in perspective, it exceeds the enthalpy of the 
cavities in the Perseus cluster 6 , Cygnus A 5 , and M87 7 by roughly 250 times, 15 times, and 
more than four orders of magnitude, respectively. Only the energy of the recently-discovered 
shock front in Hydra A 9 falls within an order of magnitude of this value. 

The bright elliptical region surrounding the X-ray cavities strongly resembles the hot 
cocoon of jet-powered radio source models 10,11 . In this interpretation, the relatively sharp 
edge of the cocoon lies at the location of an enveloping shock (Fig. 3). The two red "hot 
spots" in Fig. 4 show that the gas near the cavities is being heated by the shocks. By 
contrast, the gas surrounding the cavities in other systems such as Hydra A 4 and Perseus 6 
is relatively cool, suggesting that buoyancy, not excess pressure, is driving their outward 
advance. 

The shock properties were determined using a spherical hydrodynamic model of a point 
explosion at the center of an initially isothermal, hydrostatic atmosphere (Fig. 3). The age 
and driving energy of the shock are t s ~ 1.04 x 10 s yr and E s ~ 5.7 x 10 61 erg, respectively. 
The energy is proportional to the preshock temperature, which probably exceeds 5 keV. The 
spherical model underestimates the shocked volume, tending to underestimate total energy. 
Furthermore, since the cavities occupy a large fraction of the cocoon, they appear to be 
driving the shock, which undermines to some degree our assumption of a point explosion. 
Nevertheless, we expect the energy of the outburst to be within a factor ~ 2 of the model 
estimate, and more likely to exceed it. 

The shock energy is reassuringly close to the cavity enthalpy, and we adopt the shock 
energy as the probable value. The average jet power of the outburst is then P s = E s /t s ~ 
1.7 x 10 46 erg s _1 , comparable to a powerful quasar radiating at the Eddington limit of 
a~2x 10 8 M black hole. The flux density of the radio source at 1.4 GHz is 21 mJy 12 , 
corresponding to a monochromatic luminosity of 4 x 10 40 erg s -1 , several orders of magnitude 
fainter than powerful quasars 13 . The ratio of average jet power to monochromatic radio power 
is ~ 10 5 , enormously larger than the generally accepted factor of 10 to 100 13 ' 14 for powerful 
radio sources. Evidently, even relatively faint radio sources can be mechanically powerful 5 . 
Bright radiation from an active quasar is surprisingly absent both in the optical and X-ray 



bands. However, the cD harbors a ~ 10 42 erg sec -1 optical emission nebula extending over 
its inner 20 kpc 15 , which are commonly found in cooling flows. 

This outburst released enough energy to quench a 200 M yr" 1 cooling flow for several 
Gyr, assuming all of the energy is deposited within the ~ 50 kpc cooling region of the 
cluster. The central cooling time of the gas is roughly 1 Gyr, so the cooling that ensues 
could establish a feedback cycle of heating and cooling driven by accretion onto the central 
black hole 8 ' 16 ' 17 . A similar process occurring in other clusters would in principle maintain 
the observed levels of hot gas with short cooling times, molecular gas 18 , and star formation 19 
in cD galaxies, while preventing the development of a more massive, steady cooling flow 1 ' 20 . 
The existence of bright nebular emission located in the the cD galaxy is consistent with this 
picture 15 . 

Much of the energy is, however, leaving the cooling region, bound for the cluster's out- 
skirts. The gas mass 5.5 ± 0.7 x 10 13 M within one Mpc is being heated at the level of 
~ 1/3 keV per particle. This one event alone then provides a substantial fraction of the 
1 — 3 keV per particle of heat required to raise the entropy above the level of gravity alone 
(preheating) 21 . For a bolometric, unabsorbed X-ray luminosity of 1.1 x 10 45 erg s _1 and a 
mean temperature of 4.5 keV, the cluster departs upward in luminosity by several times from 
the relatively tight correlation between X-ray luminosity and gas temperature 22 ' 20 . Since the 
shock power exceeds the cluster's X-ray luminosity by 15 times, the shock is surely capable 
of causing this departure. The time required to radiate away the shock energy > E s /L x ~ 2 
Gyr is a substantial fraction of the age of the cluster. Therefore, this outburst will leave 
a persistent mark on the temperature and luminosity of the cluster, even after the cavities 
have disappeared. Assuming this event is not unique, substantial heating must have occurred 
recently in clusters, not just during an early preheating epoch 21 . Events of this nature would 
complicate the use of X-ray temperature and luminosity functions to probe the large scale 
structure and cosmology 22 . 

The titanic proportion of this event suggests it was powered by accretion onto a black 
hole. Equating the shock energy to O.lMc 2 gives an estimate of the minimum accreted mass 
required to power the burst of M ~ 3 x 10 8 M , itself the mass of a supermassive black 
hole. The relationship between galactic bulge luminosity and black hole mass 23 predicts a 
~ 10 9 M black hole resides there (assuming the cD's absolute visual magnitude within a 
35 kpc diameter is — 22. 4 24 ). Evidently the central black hole accreted a substantial fraction 
of its own mass in only 10 8 yr, a remarkable growth rate for such a large black hole. While 
a similar line of reasoning would apply to quasars, their ages and average jet powers have 
not been measured directly. Such a rapid rate of growth may be difficult to reconcile with 
the small scatter in the black hole mass to bulge mass relation 25 , and is at variance with the 
view that the most massive black holes have evolved slowly in the recent past 26 . 

Finally, the magnetic field strengths in clusters are typically a few /iG 27 , and evidence 
is growing for the existence of large-scale, intergalactic fields 28 ' 29 . These fields could be 
generated and dispersed by outflows from supermassive black holes 28 ' 29 . The equivalent 
magnetic field strength corresponding to the energy density within the cavities is ~ 100 /xG, 
considerably larger than the accumulated field strengths in clusters 27 . Therefore, a plausibly 



small fraction of the total energy of this one radio outburst alone channeled into magnetic 
field would magnetize the cluster. 
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Figure 1 The relationships between the X-ray, radio, and optical emission of the 
cluster. The smoothed X-ray image (left) and optical image (right) are superposed 
with the 1.4 GHz radio contours. The 40 ksec X-ray image was obtained with the 
Chandra X-ray Observatory on 1 December 2003. Approximately 75000 useful X- 
ray photons were detected. The X-ray surface brightness depressions (cavities) are 
between 10% — 20% fainter than the surrounding X-ray emission to the north and 
south of the center. Most of the cluster's emission emerges from an elliptical structure 
bounded by a shock front. We assume a flat cosmology with Ho = 70 km s -1 Mpc -1 
and Q = 0.3, corresponding to a ratio of linear to angular size of 3.5 kpc arcsec -1 at 
the redshift of the cluster throughout this paper. 

The ~ 4 arcsec resolution radio map was made with the Very Large Array telescope in 
the C configuration. The cavities are filled with radio emission. Assuming spherical 
cavities whose edges lie at the midpoints of the rims, each is roughly an arcmin in 
diameter (200 kpc) centered approximately 125 kpc to the north-east and to the south- 
west of the cluster center RA = 07 41 44.0, Dec = +74 14 38.3, (J2000). The radio con- 
tour levels are 2 x 10~ 4 x (-1, 1, 1.4, 2, 2.8, 4, 5.7, 8, 11, 16, 22, 32) Jy/beam. The radio 



contour levels on the visual image are 2 x 10 



•1,1,1.4,2,2.8,4,5.7,8,11,16,22) 



Jy/beam. The cD's i?-band surface brightness peaks at its center at roughly 19.5 
mag arcsec -2 and diminishes with radius in r 1//4 -law fashion, reaching 25.5 mag 105 
kpc). Each image is 250 x 250 arcsec (875 x 875 kpc) on a side. 
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Figure 2 Projected temperature and pressure profiles of the cluster gas. The vertical 
error bars are 90% confidence intervals, and the horizontal bars represent the bin 
sizes. The spectra used to construct these profiles were extracted from the X-ray 
image in circular apertures centered on the cD. The spectra were modeled as thermal 
emission from a single temperature plasma of uniform metallicity, attenuated by the 
foreground column of neutral hydrogen in our Galaxy. The average metallicity of 
the gas was found to be 0.4 times the solar value. The temperature and pressure 
profiles are complex. The coolest gas located in the center of the cluster is roughly 
3 keV. The temperature rises with increasing radius reaching an average of about 7 
keV between 50 and 80 arcsec. Beyond 80 arcsec the temperature drops to roughly 5 
keV, although the level of this drop is uncertain. The pressure is highest in the center 
3 x 10~ 10 erg cm -3 , and declines smoothly with radius in roughly power law fashion 
until reaching a radius of roughly 70 arcsec. The pressure there rises abruptly at the 
shock front. 
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Figure 3 Projected 0.5-7.5 keV radial surface brightness profile of the cocoon region 
compared to shock model predictions. The profile was measured in 30° sectors to the 
east and west, along the minor axis of the ellipse (PA's 90° - 120° and 270° - 300°), 
where the shock is relatively uniform. The feature at a radius of 69 arcsec (240 kpc) 
is consistent with being a weak shock. The surface brightness beyond is well fitted 
by the power law, r _/3 with (3 = 2.24 ± 0.29 (90%). To be consistent with the surface 
brightness profile beyond the shock, the gas density initially has p{r) oc r~ v , with 
i] = 1.62, and the gravitational field was chosen to make the undisturbed atmosphere 
hydrostatic. The model assumes the temperature of the unshocked gas is 5 keV. 
The Chandra 0.5 - 7.5 keV response was computed using XSPEC and an absorbed 
mekal model with foreground column density 3.49 x 10 20 cm -2 , redshift 0.216 and 
abundance 0.4 times solar (results are insensitive to these parameters). Model surface 
brightness profiles are scaled to match the observed profile in the unshocked region. 
The three lines represent model profiles for shock Mach numbers of 1.41 and 1.41±0.07 
(increasing Mach number from bottom to top). The greatest source of uncertainty in 
the age is the preshock temperature, since the Mach number is not strongly model- 
dependent (in particular, if the cocoon is axially symmetric, the age estimate is not 
sensitive to projection effects). The temperature increase for a shock of this magnitude 
is expected to be roughly 30% above the preshock value. This jump is consistent with 
the data, but the uncertainty in the pre- and post-shock temperatures are too large 
to further constrain the shock properties. The vertical error bars are 90% confidence 
intervals, and the horizontal bars represent the bin sizes. 
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Figure 4 Temperature map of the central 200 x 200 arcsec of the cluster. Redder 
colors indicate hotter temperatures. The logarithmically spaced contours, ranging in 
surface brightness from la = 1.0~ 4 cnt s _1 cm -2 pix~ 2 to 100a, show the locations of 
the cavities. The hottest regions are at the tips of the cavities, where the shocks are 
strongest. 



